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1. INTRODUCTION 


The J-integral has received much attention among researchers since it 
was introduced to fracture mechanics. As a result, a wide body of literature 
exists in this connection. The J-integral is a measure of severity of the 
deformation field at the crack tip and has proven to be a viable parameter 
for prediction of ductile crack initiation and growth under monotonic loading 
in the nonlinear regime. Because it can be determined from the far-field 
numerical solution, which is generally more accurate than the near-field 
solution, the J-integral has been quite popular in computational fracture 
mechanics. Unfortunately, however, the theoretical foundation upon which 
the J-integral is based does not permit it to be extended to more practical 
engineering problems. The path independence of the J-integral is valid only 
within the deformation theory of plasticity; hence, it cannot be used in the 
presence of substantially nonproportional loading and unloading after plastic 
deformation. The J-integral cannot be used in the presence of a temperature 
gradient and material inhomogeneity. Other factors such as body force and 
crack surface loading, although not as frequently encountered in applications 
as the above, are also excluded in the J-integral. 

In recent years, there have been considerable efforts to modify or 

reformulate the path- independent (P-I) integral such that these limitations 

•kit 

are removed. Accordingly, there emerged a number of new P-I integrals in 
the literature. Some of them are merely slight modifications of the original 
J-integral, but others are formulated on different theoretical bases. None of 
the new P-I integrals have yet been used enough to establish a consensus on 
the utility in fracture mechanics. 


* The crack surface traction can be included in the J-integral simply by 
adding the crack surface path to the integration path. 

** In the classical sense, the path-independent integral is a line integral 
performed along an open or closed curve. Here, the notion of "path inde- 
pendence" is somewhat different. As we will see, most path-independent 
integrals contain not only line integrals but also area integrals. 
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The purpose of this study is to review the new P-I integrals from the 
standpoint of theoretical basis of formulation. We will restrict our atten- 
tion to rate-independent materials and will not consider the effects of 

•fc 

inertia, body force, and large deformation. The P-I integrals in this 
category and to be reviewed in the text are the J-integral (Reference 1), 
Wilson and Yu's integral (Reference 2), Gurtin's integral (Reference 3), the 
Jg-integral by Ainsworth et al. (Reference 4), the J*-integral by Blackburn 
(Reference 5), the 5-integral by Kishimoto et al. (Reference 6), and the AT^ - 
and AT*-integrals by Atluri et al. (Reference 7). These P-I integrals will be 
recast in consistent notation, and the similarities, differences, salient 
features, and limitations will be examined. Comments will also be made with 
regard to the physical meaning, the possibility of experimental measurement, 
and the ease of computation. 


* The body force can be included easily in any P-I integral. Usually it 
generates an area integral term. 
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2. NOTATION 


In most sources reviewed in this study, the index notation has been used. 
Accordingly, we will use the index notation for all the P-I integrals. The 
variables to be used in the text are described below. Unless stated other- 
wise, the definitions of variables are as given here. For the notation 
related to integration paths and areas, the reader is also referred to 
Figure 1. 
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r r with radius e and with origin at the crack tip; 

£ e->0 is implied without limit notation 

6.. Kronecker delta 

ij 

A Increment of the subsequent variable from time T to t + Ax 

AT Atluri's P-I integral 

ATp Atluri's et al. P-I integral 

AT^ Atluri's et al. P-I integral 

e. . Strain tenor 

ij 

0 Relative temperature 

|J Lame's constants 

v Poisson's ratio 

CT. . Stress tensor 

ij 

X Time 

<]) Global thermodynamic potential 

He Imho lz free energy 
Q Specific internal energy 

(•) Elastic component of (•), such as u^, , W 

(•)^ Plastic component of (•) 

(•) e Thermal component of (•), such as e® . 



Figure 1. Integration Paths and Areas. 



3. THE J- INTEGRAL 


The introduction of the J-integral to fracture mechanics is due to Rice 
(Reference 1), although some related earlier studies can be found in Eshelby 
(Reference 8), Sanders (Reference 9), and Cherepanov (Reference 10). Later 
researchers extended the J-integral to large elastic deformation (References 
11-13). In the following, we will review the concept of the J-integral. 

Let us consider a cracked body as shown in Figure 1. We assume that the 
crack is straight and oriented in the xi-direction. We will also assume 
that the body force is negligible and the crack surface is traction-free. 
Then, the following integral denoted by J is independent of the path: 

J = J r (njW - t i u i x ) ds (1) 

where the strain energy density function W is given by 


W = 




( 2 ) 


In view that W is a uniquely defined function of the current strain, the 
J-integral is considered to be valid for nonlinear elastic materials. 


The path independence of J is easily proved by use of the equilibrium 
equation 


a. . . 




= 0 


and the divergence theorem 


(3) 


J c n i f ds = X A dA (*) 

where f is a piecewise continuously differentiable function, C is a closed 
contour, and A is the area surrounded by C. The homogeneity of the material, 
at least in the Xj direction, is also assumed to prove the path independence 
of J. In this connection, we note the work by Smelser and Gurtin (Reference 
14) who showed that the J-integral can be extended to bimaterials without 
change, provided the bond line is straight and parallel to the crack. 
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The J-integral is a measure of the severity of crack tip deformation 
when T is taken sufficiently close to the crack tip. For linearly elastic 
materials, J can be related to the stress intensity factor by direct substi- 
tution of the singular solution (Reference 15). For instance, the Mode I 
solution yields 

J = | K2 (5) 

where K = 1 for plane stress problems, and 
K = 1 - v 2 for plane strain problems 

Equation 5 has been widely used to determine from J computed from 
far-field data. 

The J-integral is also valid for elastic-plastic materials within the 
deformation theory of plasticity. For power-law hardening materials, the 
crack tip field is characterized by the HRR (Hutchinson, Rice, and Rosengren) 
singularity (References 16 and 17). The J-integral is again related to the 
strength of the singular field by 

1 

a ij = a “(^v ) n+1 V e ’ n) (6) 


and 


e. . 

ij 


- 2a 
E 


( 


EJ 

ag! n r 



e. . (0,n) 
ij 


(7) 


where ct 0 is the yield stress, n is the hardening exponent, r and 0 are the 

polar coordinates with the origin at the crack tip, 1^ is a constant that is 

a function of n (References 16 and 17), and a., and e. . are dimensionless 

ij iJ 

functions of 0 and n. A study by Kumar, et al. (Reference 18) revealed that 
J is a viable parameter for prediction of ductile crack extension under 
monotonic loading. 
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The J-integral computed from the numerical results of finite-element 
analyses, based on incremental theory of plasticity, is essentially path 
independent provided that the degree of nonproportionality in loading is not 
severe. It is also worth noting that McMeeking's work (Reference 19) shows 
the path dependence of J in the very vicinity of the crack tip where the load- 
ing is significantly nonproportional. One of the limitations on the utility 
of the J-integral is that it cannot be used in circumstances where unloading 
takes place after plastic deformation; this is significant since most crack- 
propagation problems involve unloading at either global or local level. Some 
researchers (References 20 and 21) correlated crack growth data to AJ defined 
operationally using the approximation formula for J in References 22 and 23. 
This type of approach, despite a lack of theoretical rigor, may be worthwhile 
insofar as no other parameters are available for crack-growth prediction in 
the elastic-plastic regime. 

Rice (Reference 24) has shown that the J-integral can be interpreted as 
the rate of potential energy decrease for two bodies differing in crack length 
by an infinitesimal amount, namely 

J = - dP/da (8) 

where 

P = X v W( £ij ) dV - J s t.u. dA (9) 

is the potential energy, V is the volume of the body, and is the boundary 
where traction is prescribed. In view of Equation 8, J is often referred to 
as the "crack driving force" or "energy release rate." This interpretation 
is, of course, valid within the framework of nonlinear elasticity. The rela- 
tion of Equation 8 makes possible the experimental measurement of J as dis- 
cussed by Dowling and Begley (Reference 25) and Rice, et al. (Reference 26). 
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4. WILSON AND YU'S INTEGRAL 


Wilson and Yu (Reference 2) modified the J-integral to include thermal 
strain. Their integral (for convenience, we will use J^) is given by 

J w = J r ( ni W - t.u.^) ds - ct(3A + 2 m) S a teCSe..),! -£^ 6 , 1 ] dA (10) 

where 


W 


h a. .s. . 

xj ij 


( 11 ) 


and £.. includes the thermal strain. The stress component a., is related 
xj ij 

to e. . by 


a . . = As, , 6 . . + 2 ue. . - a(3A + 2 u) 6 . .0 ( 12 ) 

xj kk ij r xj r xj 


Notice that the modification was made only for the case of homogeneous , 
isotropic, linearly elastic materials. Alternative forms of were also 
suggested in Reference 2 using different definitions of W: 


J W = J"r ^ niW * " t i u i,l^ ds + + 2fJ ) S a e 1 i 0 >i dA 


and 


J y = (n^ - t.u ) ds - a(3A + 2|j) f A 0s ii>1 dA 


( 13 ) 


(14) 


where 


W* = k o. .£. . - or(3A + 2 m)0s. . (15) 

ij ij ix 

and 

W = \ a. .£. . + a(3A + 2jj)0s . ( 16 ) 

ij ij ix 
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The path T must be replaced by T + I" c if the crack surface is not free of 
traction. The near-field expression for can be written as 

J w = Jp ( n i w " t ± u i x ) ds (17) 

£ 

Using the singular solutions (References 15 and 27), can be related to the 
stress intensity factors for a bounded crack-tip temperature. 

It is noted that material inhomogeneity can be included in in a manner 
similar to Atluri's (Reference 7); that is, the additional term in 3W/3x! due 
to material inhomogeneity can be expressed in terms of stress, strain, temp- 
erature, and the derivatives of these parameters. The resulting equation is, 
however, identical to the J*-integral (Equation 31) to be discussed later. 

The singularity of the temperature field at the crack tip, if not bounded, 
must be weaker than the square-root singularity in order for the area integral 
to exist. The possibility of an unbounded temperature at the crack tip was 
discussed in Nguyen (Reference 28) and Bui (Reference 29). Notice also that 
for a bounded 8 the integrand of the area integral in Equations 10 and 14 
behaves like r -3/2 near the crack tip while the integrand in Equation 13 
behaves like r -1/2 . In this sense, Equation 13 would be somewhat better for 
the computational purpose. 

Wilson and Yu's integral was reconsidered by McCartney (Reference 30) 
with a thermodynamic background. He showed that 

J w = Jr ^ nipv “ fc i u i,i^ ds + Ja pn6>1 dA (18 - ) 

where p is the density of the material, 4* is the helmolz free energy, r| is 
the entropy, and 0 is the absolute temperature. For the adiabatic process 
(r| = constant) , Equation 18 reduces to 

J W = Jr (aipn " t i u i,l^ ds (19) 

where Q is the specific internal energy. It is worth noting that Equation 18 
can be written as 
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J w = h ’ t i U i,l ) ds " *^A p6n>1 dA 


( 20 ) 


Here, 4* = Q - 0r) has been used. Then, replacing T with T + F c and using the 
divergence theorem, the near-field expression is obtained as 


J w = Jr (nipQ ‘ t i U i,l ) ds 

c 


( 21 ) 


It is of interest to notice the work of Nguyen (Reference 28) and Germain, et 
al. (Reference 31) in which they showed that Equation 12 can be expressed as 


J w = - 34>/9a (22) 

where <|> is a global thermodynamic potential with arguments such as strain, 
temperature, crack length, and other internal parameters. Equation 22 is 
equivalent to Equation 8 for J and, as such, carries a significant meaning 
from the standpoint of experimental measurement. Notice also that Equation 2 
is equivalent to Equation 17 if fi is identical to W at the crack tip. This 
condition is satisfied if the entropy at the crack* tip is maintained constant. 
For a stationary crack, this implies a finite crack tip temperature. 
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5. GURTIN'S INTEGRAL 


Gurtin (Reference 3) proved the following conservation law for the 
thermoelastic-field-satisfying Equations 3 and 12 and the steady-state 
temperature condition, 0,.. = 0: 


n.W 

x 


t k U k , i 


« . 2 (3M 2tj). 2 . 0 2 n + 

2(X+JJ) ° “i (X+M) 



30 \ 
U i 3n ) 


dA 


0 (23) 


where S is the boundary of the volume, 3/3n = n. 3/3x., and W is the strain 

J J 

energy density defined by 


W = " £ ij £ ij + 2 (£ tk )! 


(24) 


Equation 23 is valid in two or three dimensions. In particular, for the 
former with traction-free crack surface, 



n]W - t k u 


k,l 


a 2 (3\+2y) 2 n , ay(3X+2u) /. 3ui 
2 (\+m) 9 ni + (X+m) [ 9 dT 


ui 


30 \ 
3n / 


ds 


(25) 


is path independent provided that 0ui, 2 and Ui0, 2 are continuous across the 
crack. This condition is a severe restriction from the application stand- 
point. For instance, the continuity of the two quantities cannot be satis- 
fied for Mode I crack problems. It would be more practical to change the 
integration path from T to T + T , thereby eliminating the restrictions. 

With this modification and with the assumption that the singularity of the 
temperature at the crack tip is weaker than r -1/2 , the near-field expression 
for J„ can be written as 

Jq = J r (n x w - ti^ 1 ) ds (26) 


This equation and Equation 17 yield the same results, although different 
definitions of W were used. 

Equation 25 clearly has a computational advantage over because it does 

not contain area integrals, and it is possible to determine experimentally 
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by measuring the variables in the integrand along the boundary. This type of 
measurement was performed by Read and McHenry (Reference 32) to determine J 
for a single-edge notch specimen. Notice, however, that the J^-integral does 
not require the steady-state temperature condition, but does. 


% 
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6. THE J Q - INTEGRAL BY AINSWORTH, ET AL. 


The Jg-integral introduced by Ainsworth, et al. (Reference 4) is given 


by 


J Q = Jr ( n i w " t.U. ,) ds + /. CT. .£®. - 

0 J r v 1 11,1 j a i] ij,i 


dA 


(27) 


where e®. is the thermal strain, 
ij 


£ . 

W(e! .) = J 1J a. . de'. 

V IT 7 J 'll 1 ’ 


IJ 


IJ IJ 


(28) 


and 


£’. . 

ij 




Notice that e! . includes the elastic and plastic strain due to mechanical 
ij 

loads. For elastic deformation it can be easily shown that Jg is identical 
to Jy The J-integral is valid only within the deformation theory of plas- 
ticity, and it cannot be used with unloading after plastic deformation and 
material inhomogeneity. Ainsworth, et al. (Reference 4) asserts that Jg is 
related to the potential energy by 


J 


e 


= lim r— [/_ n .CT. .Au. ds - f, CT. .Ae. . dA] 
Aa+0 Aa r J 1J 1 A 


(29) 


where A denotes the increment as the crack length changes from a to a + Aa. 
In fact, they started from Equation 29 and arrived at Equation 27 in their 
derivation. It is noted, however, that Jg given by Equation A7 in Reference 
4 is not identical to Equation 27 unless the integration of njW along T 

A ^ 

vanishes. Notice also that their Equation A7 is the same as the J-integral 
given by Equation 36. 
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7. THE J*-INTEGRAI BY BLACKBURN 


Blackburn (Reference 5) proposed the J*-integral defined by 

J* = fr- (%o. .u. . dxo - t.u. , ds) 

J r c v ^ij i,j 2 n,l 


(30) 


or alternatively 


J* = JV,,- (IfcCT. .u. . dx 2 - t.u. , ds) + f. (^a. .u. ..-^ct.. _u. .) da (31) 
J r+r ij 1,J 2 1 1,1 J A v ^ij i,jl 1J,1 i,J 


It can be shown that the second term in Equation 31 vanishes in the elastic 
case; thus, J* is identical to J. In the thermoelastic case, it is identical 
to J^, Jg, and Jg in view of the near-field expression of these P-I integrals. 
Furthermore, the J*-integral given by Equation 28 is path independent even in 
the presence of material inhomogeneity since only the equilibrium and the 
divergence theorem were used to obtain Equation 31 from Equation 30. For 
power-law hardening materials, direct substitution of the HRR field furnishes 
the relation between J* and J: 


J* = g(n)J (32) 

where n is the hardening coefficient, and g(n) can be computed once the 
angular variation of field variables is known (see Appendix). For general 
elastic-plastic materials with smooth stress-strain curves, the integrands in 
Equation 31 are continuous everywhere in the body if the material undergoes 
either loading or unloading in the global sense. For loading conditions that 
give rise to loading and unloading zones in the body at the same time and for 
materials with kinks in the stress-strain curves (such as bilinear material), 
the continuity of the stresses and the derivatives of displacement must be 
guaranteed across the loading-unloading boundary or the elastic-plastic 
boundary for the path independence of J*. The continuity of the stresses can 
be proven from the equilibrium of a thin slice of the body containing the 
boundary and from the continuity of the effective stress. It is also believed 
that the strain will be continuous across the boundary (excluding the elastic/ 
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perfectly plastic material) considering that the stress continuity will ensure 
the continuity of the elastic strain and the plastic strain will change con- 
tinuously across the boundary (but the plastic strain rate is discontinuous). 
Continuity of strain will then guarantee the continuity of the displacement 
gradient; hence, use of the divergence theorem (Equation 4), and consequently 
the path independence of J*, is also justified under these circumstances. Un- 
fortunately, no numerical results for J* are found for problems of this sort. 


The major difference of J* from other P-I integrals is that the strain 
energy density function does not appear in the integral. Instead, an explicit 
expression SsCT^u^ . was used. This type of approach will expand the scope of 
the applicability of the integral, but it makes the physical meaning obscure. 
In fact, the physical meaning of J* in the plastic regime is not known and may 
be difficult to ascertain. The determination of J* through experiment is not 
likely unless J* is related to potential in a manner analogous to Equation 8. 
Notice that Equation 30 requires accurate measurement of near-field variables, 
and Equation 31 has an area integral that can be determined only by measuring 
the field data everywhere in A. Neither of these is expected to be feasible. 
Other P-I integrals containing area integrals have the same problem. 


The computation of Equation 31 requires caution in view that the inte- 
grand behaves like r -2 in the vicinity of the crack tip. The second term can 
be written near the crack tip as follows: 


r 11 

J J" [r -2 f (0 ) + less singular terms] r d0 dr (33) 

0 “7T 


Thus, for J* to exist (this is guaranteed by Equation 27), we must have 


/ 


-it 


f (0) d0 = 0 


(34) 


This can be readily verified for linearly elastic materials, as noted 
by Atluri (Reference 33). For other types of materials, the angular 
distribution of the field variables is not known explicitly. 

Blackburn, et al. (Reference 34) presented some numerical results of 
J* for specimens subjected to loading, unloading, and temperature gradient. 
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8. THE J- INTEGRAL BY KISHIMOTO, ET AL. 


The J- integral proposed by Kishimoto, et al. (Reference 6) is given by 


J = -J r t.u. ds (35) 

end 1 1,1 


where r en( j is the path surrounding the so-called "fracture process zone" m 
which the continuum mechanics fails to work effectively. In the process to 
arrive at Equation 35, the equilibrium equation and the divergence theorem 
were used. Also, it was assumed that the displacement field on I~ en( j is main- 
tained constant in the crack propagation, namely 3u^/3a = 0 on With 

this assumption, 3 represents the rate of work done on the fracture-process 
zone by the surrounding medium. Then they convert Equation 35 into the 
following by assuming that r en( j shrinks and vanishes: 




r . _ t.u. , 

r+r c i 1,1 


ds + f. ct. . e. . , dA 
J A ij ij,l 


They further assumed that 


(36) 


f r W e ds = 0 (37) 

e 

and obtained 


a , 

J = J, 


(niW - t.u. .) ds + J. a. .£*. . 
r+r c v 1 i i,l' J A ij ij , 1 


dA 


(38) 


where 


e*. 

ij 





+ e«. 

ij 


( 39 ) 


and 
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( 40 ) 



a. .ds?. 
ij 


For linear elastic problems, Equation 37 is not satisfied unless a non- 
singular crack-tip field is enforced in the analysis. With the singular 
elastic field. Equation 38 is equivalent to 


3 = / r (nxW 6 - t i u £ j) ds (41) 

£ 

which is the usual J-integral. Thus, 3 computed from Equation 38 is not 
identical to 3 computed from Equation 36. For thermoelastic problems it can 
be easily shown that 3 given by Equation 38 is identical to Jg, J^, Jg, and J* 
provided that the material is homogeneous. It is also noted that the physical 
meaning they claim is valid if there is a finite, rigid, fracture-process zone 
at the crack tip - independent of the crack size. It is unclear, however, 
whether or not the singular crack tip field is recovered when r ^ shrinks and 
vanishes and, if it is recovered, whether or not the physical meaning remains 
intact. 

For elastic-plastic problems, the crack tip field is dominated by plastic 
deformation, as verified by Hutchinson (Reference 16) for power-law hardening 
materials. The hypothesis of Equation 37 is, therefore, acceptable for a 
vanishingly small fracture-process zone, and Equation 38 would be practically 
identical to Equation 36. It is also mentioned that 3 , given by Equation 36, 
can be related to J for linear elastic and power-law hardening materials by 


3 = 


/ J(2v 2 + v + 3)/4 
J(3 - 2v)/4(l - v) 
J h(n) 


elastic materials, plane stress 
elastic materials, plane strain (42) 

power-law hardening materials 


where h(n) is a function of the hardening coefficient only (see Appendix ) . 

The J-integral given by Equation 36 can be used for elastic-plastic 
materials subjected to loading and unloading. Aoki, et al. (Reference 35) 
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demonstrated the path independence of J for a specimen subjected to global 
loading and then unloading. For problems undergoing local unloading, no 

A 

numerical results are available; however, it is expected that J will again be 
path independent in view that the integrand of the line integral will be con- 
tinuous across the loading/unloading boundary as described in the discussion 
of the J* integral. Equation 36 can also be used with a temperature gradient 
and material inhomogeneity since no constitutive relations were assumed in the 
derivation. 

a , 

The computation of the area integral in J is similar to that of the J* 

integral. The experimental measurement of 5 would be as difficult as that of 

J*. 
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9. THE AT- INTEGRALS BY ATLURI, ET AL. 


Atluri (Reference 36) proposed a path-independent integral, given in the 
incremental form, which takes into account finite strain, body force, inertia, 
and crack surface traction. Temperature gradient and material inhomogeneity 
were not included in the formulation. The small-deformation version of the 
AT-integral is given by 

“ = • l r>r c I “ lAM ‘ a j (o jk + Ao jt )A \,i 1 ds - 4 0 jk,i il, j,k dA (43) 

= R + JV [niAW - n.(a.. + Act..)Au. .] ds 
r_ l 1 j v ij ij' i,l J 


where 


AW = (a. . + %Aa..)Au. . 
ij ij i»J 


(44) 


and R is the contribution of the discontinuities in the material response 
along the loading/unloading boundary. The incremental stress AC4 ^ is related 
to the incremental strain by the Prandtl-Reuss equation: 


12a|J 2 AE o' ff! 

Act.. = 2pAe . . + \S .Ae. . ■ ■ - ■ (45) 

1J 1J 1J kk CT' CT' (6\) + 2dI/d\P) 

mn mn 


where 


ct'. . 
ij 



1/3 6 . . ct 

ij 


kk 


F = l ct’ . CT’. . 2 - F 0 (W P ) = 0 

2 ij ij ° 


and 


(46) 

(47) 


P 

W P = / « 


. P 
ct. . de. . 
ij ij 


( 48 ) 
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The AT-integral defined by Equation 43 includes the measure of the severity 
of conditions near the crack tip and also the effect of the transition from 
loading to unloading zones in A. The AT-integral was proven to be path inde- 
pendent for paths in either the entirely loading zone or entirely unloading 
zone. Atluri (Reference 36) also showed, for the case of nonsingular crack 
tip field, that the AT-integral implies the rate of change in potential energy 
(per unit crack length) in the time interval x to T + At of two bodies that 
are identical in shape and load history but have an infinitesimal difference 
in crack length. 


Atluri, Nishioka, and Nakagaki (Reference 7) later presented two modified 
versions of the AT-integral. In the new integrals, denoted by AT^ and AT^, 
the contribution of the elastic-plastic boundary was expressed in terms of 
field variables. The AT^ and AT^ integrals are given by: 


AT* = ; r [ Qi aw 


(t. + At.)Au. , - At.u. .] ds 

l l 1,1 l 1,1 



[n.AW 

l 


(t. + At.)Au. - 

l l 1,1 


At.u. ,] ds 
i i,l J 


(49) 


+ f. [Act. . (e . . - + . . , ) - Ae . . (a. . , + \ha. . . ) ] dA 

A ij ij,l ij,l ij ij,l ij,l 


and 


AT = f 0 [niAW - (t. + At.)Au. . - At u. ,] ds 
p J S 1 l i ' i,l l i,l J 


= J r+r [ ni AW - (t + At.)Au - At u ] ds 
c ’ ’ 


(50) 




dA 


In the AT -integral, S is the external boundary including the crack surface, 

Aj, is the area enclosed by T, and A g is the total area. The area integrals 
on the extreme right-hand sides of Equations 49 and 50 vanish when the loading 
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is proportional. Thus, inasmuch as the deformation theory of plasticity is 
applicable, the following identity holds: 


AT* = at = AJ 
P P 

where 


(51) 


AJ = J r [njAW - (t., + At i )Au i l 

C ^ 


At.u. .] 
1 i,l J 


ds 


(52) 


is considered the increment of the J-integral due to the increment of the 
external loading in the time interval X to T + At. 

Equation 49 shows that AT^ is a direct measure of severity of deformation 
at the crack tip, but it is not experimentally measurable. On the other hand, 
in view of Equation 50, AT^ is not related to the immediate crack-tip field 
unless the loading is proportional. Notice, however, AT^ can be determined by 
measuring field data along the external boundary. Nakagaki, et al. (Reference 
37) presented some numerical data to demonstrate the path independence of AT^ 
and ATp for a compact tension specimen that was plastically loaded, unloaded, 
and then reloaded. They also computed AJ at the peaks using Equation 52 and 
showed that it is no longer path independent after the specimen is unloaded. 

A salient feature of the Al* and AT^-integrals is that they are based on 
the incremental theory of plasticity; this enables these P-I integrals to be 
used even in nonproportional loading and elastic unloading following plastic 
deformation. It can also be verified, by adding -(2p + 3\) aa^AQ to Equation 
45 and following the procedure in Reference 7, that Equations 49 and 50 can be 
used without change for thermomechanical loading. Material inhomogeneity is 
also taken into account in Equations 49 and 50. 

The computation of AT* would be similar to other P-I integrals involving 
area integrals. The behavior of the integrand of the area integral of AT* in 
the vicinity of the crack tip is similar to the area integral of Equation 31. 

Atluri et al. (Reference 7) showed that AT is related to the rate of an 

P 

incremental energy (All) with respect to crack length if the loading is propor- 
tional. For nonproportional loading, the physical meaning of AT^ is not clear 

nor is that of AT*. 

P 
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10. ENGINEERING APPLICATIONS 


In this report, several P-I integrals have been examined on a theoretical 
basis. To be useful for engineering applications, a P-I integral must satisfy 
each of three conditions. 

First, the integral must be path independent for realistic loading and 
temperature conditions. The integral must be calculable without difficulties 
and must be reasonably path-independent; that is, path- independence must be 
achieved at least numerically, even though no theoretical justification 
exists, outside of the intense deformation zone at the crack tip. 

Second, the integral must be determined (or closely approximated) from 
measurements made on laboratory test specimens at various stages of loading. 

If the integral can be expressed as the rate of change of a potential with 
crack length, it can in principle be determined from a load-displacement 
record. If the integral consists only of a line integral applicable over the 
specimen boundaries, it may be measured using strain gages and extensometers 
as in Reference 25. Area integrals would require a whole-field stress and 
strain analysis. While such may be possible, they are tedious, and results 
are not immediately available. 

Third, the integral must correlate various types of crack-propagation 
behavior. That is, crack propagation must depend on the imposed value of the 
integral and be independent of crack size and specimen or component geometry. 
This condition is most likely to be satisfied if the integral is a measure of 
the crack-tip severity and if a sound physical meaning is associated with the 
P-I integral in relation to crack propagation. 

The above three aspects are summarized in Table 1 in conjunction with 
the deformation and temperature conditions under which the integrals are 
path independent. 


22 



Table 1. Summary of P-I Integrals. 


Physical 

Measure U S S i iSa^ L 

ll ' of Crack 

P-1 Tip Thermo- 

Integral Severity Elastic elastic 


Capability to Handle 


Nonprop. 

Prop. & Loading 
Plastic loading Unloading 
(plastic regime) 


Thermal 

Strain 


Material 

Inhomogenlty 


Computation 

(integrals 

involved) 


Experimental 

Measurement 


J 

J w 

Yes 

Yes 

_3P 

'aa 

jsP 

"aa 

_a*( 3 ) 

"aa 

aP 

"aa 

Yes 

No 

J G 

Yes 

_ aP 
” aa 

.11 

‘aa 

- 

No 

J e 

Yes 

aP 
" aa 

aa 

"aa 

_ aP 
" aa 

Yes 

J* 

Yes 

. 8P 
aa 

"fa Uhknown 

Yes 

5 

Yes 

Rate of work done to crack 
tip by surrounding 
material (4) 

Yes 

< 

Yes 

r-gj for prop, loading 


Yes 

aT p 

No 

Unknown 
L loading 

for non prop. 


Yes 


No 

No 

No 

Line 

Yes 

No 

Yes 

No 

Line 
♦ Area 

Yes 

No 

Yes 

No 

Line 

Yes 

No 

Yes 

No 

Line 
♦ Area 

Yes 

Yes 

Yes 

Yes 

Line 
♦ Area 

No 

Yes 

Yes 

Yes 

Line 
a Area 

No 

Yes 

Yes 

Yes 

Line 
♦ Area 

No< 6 > 

Yes 

Yes 

Yes 

Line 
♦ Area 

Yes 


Note : (1) J ■ Jy ■ J G ■ J 0 ■ J* for elastic and thermoelastic deformation of homogeneous materials, 

J ■ EAT* ■ EATp for monotonic proportional loading. 

(2) P ■ Potential energy, a ■ Global thermodynamic potential , An ■ Incremental energy 

(3) Provided that the Internal energy Is Identical to the strain energy at the crack tip 

(4) Kith the assumption of a rigid fracture process zone at the crack tip Independent of the crack size 

(5) Yes If It can be expressed as the rate of a potential, or if It can be expressed as a line Integral 
along the boundary 

(6) Yes if the loading Is proportional 
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11. CONCLUSIONS 


We have reviewed a number of P-I integrals proposed for use in elastic- 
plastic fracture mechanics. The results are summarized in Table 1. It is 
seen that there exists as yet no P-I integral with well-established physical 
background and all the desirable features encountered in applications. 

The J-integral and J^-integral may be used only for monotonic loading 
without substantial nonproportionality. For cyclic loading, utilization of 
these P-I integrals may be possible only by use of appropriate operational 
definitions of the range of these quantities. 


The Jy and Jg integrals are usable only for thermoelastic problems with 
homogeneous material properties. These integrals may be useful for prediction 
of crack growth in a rather small temperature gradient field and under small- 
scale-yielding conditions. 


The path independence of the J*, 3, AT^, and AT^ integrals is maintained 
for more general elastic/plastic problems including nonproportional loading, 
unloading, temperature gradients, and material inhomogeneities. This clearly 
is a salient feature upon application to crack-growth analyses under cyclic 
and thermomechanical loading; however, the physical meaning of these P-I 
integrals should be further investigated. 


Finally, it is noted that no significant effort has yet been made to 
utilize these P-I integrals to consolidate crack-growth data under mechanical 
or thermomechanical load cycling. Experimental efforts along with further 
analytical studies are needed to evaluate the utility of these P-I integrals 
in fracture mechanics. 
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APPENDIX 


Substitution of the HRR singular solution into the near-field expression 
of the J* integral furnishes 


J r (%a. .£. -nj - a. .n.u. ) ds = aK n+1 e^ n+1 ^ s 2 ^ +1 I* 

T_ v ij xj ij j i,l 


(53) 


where s = (2n + l)/(n +1), e is the radius of T , K is a constant related to 

w 

the boundary condition (Reference 16), and 




#wn+1 _ ** <V <*«/ A/ 

% ct cos0 - sine[CT r (u 0 - 3u r /30) - a rQ (u r + 3u 0 /30)] (54) 


+ [n(s - 2) + 1] cos0 ( CT r u r + 5 r0 "W J d8 

/V /W A/ A/ 

Here CT g , , a 0 , a r0 , u r , and u 0 are functions of 0 only (Reference 16), 

associated with the angular distribution of the corresponding stress or 
displacement components in the vicinity of the crack tip. Notice that I* 
is a function of n only. Equation 53 can be written as 

J* = g(n) J 

where g(n) = I*/I, and I is defined by Equation 24 of Reference 16. 
Similarly, 

J = h(n) J 


where 


h(n) = I/I 


and 
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I = - J (sin6[cr r (u 0 - 3^/36) - CT r0 (u r + 3u 0 /38)] 


+ [n(s - 2) + 1] cos6 (cr r u r + S rQ u 0 )^ 


d0 


( 55 ) 


Again, I is a function of n only. 
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